Incompressibility of finite fermionic systems: stable and exotic atomic nuclei 
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The incompressibility of finite fermionic systems is investigated using analytical approaches and 
microscopic models. The incompressibility of a system is directly linked to the zero-point kinetic 
energy of constituent fermions, and this is a universal feature of fermionic systems. In the case of 
atomic nuclei, this implies a constant value of the incompressibility in medium-heavy and heavy 
nuclei. The evolution of nuclear incompressibility along Sn and Pb isotopic chains is analyzed using 
global microscopic models, based on both non-relativistic and relativistic energy functionals. The 
result is an almost constant incompressibility in stable nuclei and systems not far from stability, 
and a steep decrease in nuclei with pronounced neutron excess, caused by the emergence of a soft 
monopole mode in neutron-rich nuclei. 

PACS numbers: 21.10.Re, 21.65.-f, 21.60.Jz 



I. INTRODUCTION 



Recent advances in experimental and theoretical inves- 
tigations of nuclear incompressibility have also enabled 
studies of the evolution of incompressibility in isotopic 
chains [lH6|. As for many other structure phenomena, 
for instance the evolution of magic numbers [3] , a study 
of the isotopic dependence might provide a more gen- 
eral and deeper insight into the inner workings. It has 
recently been inferred that measurement of isoscalar gi- 
ant monopole resonances (ISGMR) may probe nuclear 
incompressibility in a range of low nucleonic densities, 
rather than exactly at saturation density [8|. Experi- 
ments on isotopic chains 0, HIj possibly also extended 
to exotic nuclei [^, [ifll , are therefore important for con- 
straining the density dependence of the incompressibil- 
ity. Investigations of isotopic chains also enable studies 
of pairing and shell effects on incompressibility. It has 
been shown that one cannot neglect the effect of shell 
closure on nuclear incompressibility, and pairing correla- 
tions affect the position of the GMR excitation energy 
[J,[a|. The effect of shell closure (and therefore magicity) 
on the incompressibility has been discussed in Ref. [ll| . 
However, there is a general consensus that pairing ef- 
fects alone are not large enough to explain the available 
data [y, [T^, [IjI- Early microscopic studies of the nu- 
clear incompressibility were carried out using Skyrme ef- 
fective forces. Later, with the inclusion of relativistic 
mean field models, detailed comparisons between differ- 
ent model predictions led to the conclusion that the nu- 
clear matter incompressibility extracted from e.g. ^"^^Pb 
lies in the range « 240±20 MeV (cf. [ij] and references 
therein). The remaining uncertainty reflects our incom- 
plete understanding of the specific effects related to finite 
nuclei, and further justifies the need of studying extended 
isotopic chains. 



The theoretical framework of energy density function- 
als (EDF) [15,, ,16,] is the tool of choice for a quantitative 
analysis of these phenomena. Preliminary studies, un- 
dertaken in the 1980's, analyzed the dependence of the 
nuclear incompressibility on mass number, but only four 
doubly magic nuclei were considered and it was suggested 
that the nuclear incompressibility is rather constant for 
heavy nuclei [13, [3 • Another aspect is that the occur- 
rence of low-energy monopole strength in exotic nuclei 
has recently been predicted |19l - [2H . and this might af- 
fect the nuclear incompressibility. More generally, the 
present study of extended isotopic chains could also be 
useful for the understanding of incompressibility in other 
fermionic systems, or for astrophysical processes such as 
core collapse supernovae. The use of both analytical and 
fully microscopic approaches will provide complementary 
information about the incompressibility in both stable 
and exotic nuclei. 

Section [IT] defines the analytical relations between the 
incompressibility and other relevant nuclear quantities. 
Section IIIII compares the analytical results with the mi- 
croscopic ones obtained using the EDF framework. A 
brief summary and conclusions are given in section IIVI 



II. INCOMPRESSIBILITY OF A FINITE 
FERMIONIC SYSTEM 



In this section we review the relations between the nu- 
clear incompressibility K^i and the nucleonic kinetic en- 
ergy in the following analytical models: infinite nuclear 
matter, Fermi gas, and the harmonic oscillator. For in- 
stance, in the free Fermi gas approach the incompressibil- 
ity of the system is directly linked to the Fermi energy 



A. Definitions 

The incompressibility of a fermionic system composed 
of A particles is defined as the second derivative (curva- 
ture) of its energy with respect to the mean square radius 
at the minimum Il7t: 
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where the label g.s. indicates that the quantities are eval- 
uated in the ground state. For the isoscalar monopole 
transition operator 



i=l 

the dielectric theorem [13, [H, [H yields 
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where ?7i„i is the inverse energy weighted sum rule, pro- 
portional to the polarizability of the system. 

The energy of the isoscalar GMR can be defined by the 
following relation [13, El]: 
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The moment nii is evaluated by the double commutator 
using the Thouless theorem 23] : 
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where m is the fermion mass, that is, mc^ « 938 MeV 
for the nucleon. Eq. ©,(111) and © lead to ^: 
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B. An analytical study 

An analytical derivation of the nuclear incompressibil- 
ity K^ in the nuclear matter, the Fermi gas, and the har- 
monic oscillator models, will provide a complementary 
insight to the microscopic results of the next section. 



1. Nuclear matter and the Fermi gas 

Let us first consider infinite homogeneous nuclear mat- 
ter. The nuclear matter incompressibility at saturation 
Kao is understood fairly well, whereas the behavior of 
the, more general, density-dependent incompressibility 



has not been studied much. By defining the general 
density-dependent incompressibility K{p) [sl, |24| 
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and neglecting inter-nucleon interactions, one finds |17l . 



i^(p) = 6.^ = ^ — 
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This equation exhibits the density dependence of the 
incompressibility, known to play an important role in 
the relation between the measurement of GMRs and the 
equation of state [8|]: in the case of non- interacting nu- 
clear matter the density-dependent incompressibility is 
proportional to p^^^. At saturation density (po=0.16 
fm-3) Eq. © yields K{pq)=220 MeV, and at the 
so-called crossing density Q (pc—OA fm^'^) one gets 
K{pc)=160 MeV. This shows that the density depen- 
dence of the incompressibility of non-interacting nuclear 
matter presents a good approximation at saturation den- 
sity but not at the crossingdensity, for which Kc —40 
MeV is the expected value [g]. This can be understood 
because at the crossing density interaction effects are ex- 
pected to be more important, whereas at saturation den- 
sity the homogeneous nucleonic system can simply be 
modeled by non-interacting hard spheres. Considering 
an average spacing tq among the constituents, one gets 
for the homogeneous density of the system: 
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In the case of nucleons, ro «1.2 fm yields po ~ 0.15 fm^"^, 
and this shows the validity of the model in the vicinity 
of saturation density. 

Approximating the nucleus by a free Fermi gas of nu- 
cleons, we study a spherical fermionic A-body system 
with a typical equilibrium inter-nucleon distance tq . The 
radius of the system is then R=roA^'^. To establish a 
link between the incompressibility Ka and the zero-point 
energy, one inserts Eq. ^ in Eq. ©, and this leads to 



Ka = ^(9^)2/3To « 7To 



(10) 



where Tq is the zero-point kinetic energy (ZPE) of a 
fermion in the system [25[: 
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The ZPE, first discussed by Mottelson in Ref. 125|, 
depends only on the characteristics of the constituent 
fermions, but not on the system itself. For nucleons Tq s» 
30 MeV, and this yields Ka « 210 MeV. It should be 
noted that Blaizot also derived a relation between Ka 



and the average kinetic energy [17|. However, the aver- 
age kinetic energy depends on the model and the system 
considered. Here the incompressibihty is hnked to the 
ZPE, which does not depend on any model nor on the 
specific system. Eq. (|10p shows that the larger the ZPE, 
the system becomes more difficult to compress: quan- 
tum localization of individual nucleons determines the 
nuclear incompressibihty. As a consequence, Ka should 
be independent of A because the incompressibihty is di- 
rectly linked to the nucleonic ZPE (provided, of course, 
that the surface to volume ratio is small enough). 



2. Harmonic oscillator 

In the harmonic oscillator (HO) approximation the 
L = radial compression mode (GMR) can be repre- 
sented by a 2hbj mode [2q|, where hw is 
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Assuming spherical density, (r^) in Eq. ^ can be ap- 
proximated by 3-R^/5 where R is the nuclear radius, and 
the GMR energy becomes: 
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In the HO approximation £'GMR=^2fi,aj, and from 
Eqs. dnj and ((El) one obtains: 
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Therefore, in the HO approximation the compression 
modulus is also found to be constant and directly re- 
lated to the zero-point kinetic energy (ZPE). This again 
means that the higher the ZPE, the harder it becomes 
to compress the system. For nuclei Ka is expected to be 
a constant in the lowest-order approximation, depending 
only on the nucleon mass and the typical inter-nucleon 
distance. Eq. ([Tl)) yields Ka — 140 MeV, and one ex- 
pects that the HO value is more accurate than the one 
obtained in the Fermi gas model. This is because the HO 
potential presents a much better first-order approxima- 
tion to the nuclear mean-field potential. 

If independent of A, however, the value Ka ~ 140 MeV 
does not correspond to the limit of infinite nuclear matter 
when A — >■ cx), that is, iiroo=220 MeV. The reason is that 
the HO value has been derived under the assumption that 
the GMR is a 2^10; mode, but this is not valid any more 
for infinite nuclear matter. In that case, it is the incom- 
pressibihty of the Fermi gas that should be used, yielding 
iiryi=210 MeV as described above. This difference is also 
confirmed by the fact that: i) Ka ~ 140 MeV ispredicted 
by microscopic calculations in finite nuclei [5|, |6| , and ii) 
Koo > 200 MeV is obtained by microscopic extrapolation 
of GMR measurements to infinite nuclear matter [8i . 



III. MICROSCOPIC CALCULATION OF THE 

NUCLEAR INCOMPRESSIBILITY IN ISOTOPIC 

CHAINS 

The first fully microscopic studies of the behavior of 
K^ as a function of the number of nucleons were per- 
formed in the 80's [13, lla]- At that time only four spher- 
ical magic nuclei were considered: ^^O, ''''Ca, ^"Zr and 
^'^^Pb. It was inferred that Ka is not very sensitive to A, 
although a general conclusion was hindered by the small 
number of nuclei included in the analysis. One should 
also emphasize that light nuclei are unlikely to fit into 
the general picture because of the large surface to volume 
ratio. In light nuclei the experimental ISGMR strength 
distribution displays pronounced fragmentation, and it 
is hard to extract a well-defined value for the monopole 
energy. 

To investigate the evolution of nuclear incompressibil- 
ity along isotopic chains, it is necessary to use global, 
EDF-based models, that also take into account pairing 
effects. The constrained-HFB (CHFB) and the quasi- 
particle random-phase approximation (QRPA) [27[ are 
employed in the present study. For the latter, both rela- 
tivistic and non-relativistic implementations will be used 
to reduce possible model dependencies. Our aim is to 
verify, in a microscopic calculation, the constant behav- 
ior of Ka predicted by the analytical model described in 
the previous section. 



A. Formalism 

Microscopic EDF-based calculations of isoscalar GMRs 
are usually performed using the CHFB or the QRPA 
frameworks [27|. In this study we calculate the GMR 
excitation energies using the CHFB approach based on 
a Skyrmc EDF, and QRPA calculations are carried out 
using both Skyrme and relativistic functionals. 



1. Constrained-HFB 

The extension of the constrained HF method [l^ to 
include pairing correlations, that is, the CHFB approach 
has been detailed in Ref. [20] , and applied also in [3| . In 
the CHFB the energy of the GMR is computed using Eq. 
(|3]), with the TOi moment determined by the HFB r.ni.s. 
radius (Eq. [5]). To evaluate the 7tt._i moment, the con- 
straint associated with the isoscalar monopole operator 
is added to the HFB Hamiltonian 
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The ?7i_i moment is obtained from the derivative of the 
expectation value of the monopole operator in the CHFB 
eigenstate |A), using the dielectric theorem 
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In the present calculation we employ the Skyrme func- 
tional SLy5 [2^, together with the density-dependent 
delta-pairing interaction 
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with yo=-285 MeV fm^, 7?=0.35, a=l and po=0.16 fm^^. 
This pairing interaction has been designed to repro- 
duce the experimental two-neutron separation energies 
for the even-even isotopes ^^^"^■^*Sn [2^. The energy 
cutoff (quasiparticle energy) is Emax = 60 MeV, and 

Jmax — LofZ. 



2. Quasiparticle Random Phase Approximation 

For completeness we briefly outline the details of the 
QRPA calculations. The Skyrme HF-BCS equations for 
the ground state are solved in coordinate space. The ra- 
dial mesh extends to 18 fm, with a spacing of 0.1 fm. 
For all nuclei analyzed in the present study we have ver- 
ified the stability of the results with the variation of the 
size of the mesh. Pairing correlations are taken into ac- 
count by a contact volume-pairing interaction, with the 
strength Vb = -218(-265) MeV fm^ for Sn (Pb) iso- 
topes. The pairing strength is adjusted to reproduce the 
empirical pairing gaps extracted from data on odd-even 
mass differences (five-point formula). The single-particle 
continuum is discretized by placing nuclei in a spherical 
box of radius 18 fm. For each set of spherical single- 
particle quantum numbers {I, j), the QRPA model space 
includes unoccupied states up to the maximum number 
of nodes given by Umax — ?^iast+12, where niast is the 
number of nodes for the last occupied state with a given 
(/, j). The convergence of the results has been verified 
by calculating the strength function of the ISGMR. For 
a state with good angular momentum and parity J^ , the 
QRPA matrix equation reads 
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where fta;„ is the energy of the n-th QRPA state. X" and 
Y" are the corresponding forward and backward ampli- 
tudes, respectively. The explicit expressions for the ma- 
trices A and B are given, for instance, in Refs. |30l[31|. 
The particle-hole {p — h) matrix elements are computed 
from the Skyrme energy density functional including all 
terms, that is, also the two-body spin-orbit and two-body 
Coulomb interactions. 

The resulting moments of the QRPA strength distri- 
bution are calculated as 
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where S{E) = E„ l(0|'3l")P'^(-E^ - ^n) is the strength 
function associated with the monopole operator defined 
in Eq. @. 



In the present analysis we also employ the QRPA 
based on relativistic nuclear energy density function- 
al. More details about this framework are given in 
Refs. [32h34| . It is realized in terms of the fully self- 
consistent RQRPA based on the relativistic Hartree- 
Bogoliubov model (RHB). The RQRPA equations are 
formulated in the canonical single-nucleon basis of the 
RHB model, by employing effective Lagrangians with 
density-dependent meson-nucleon couplings [33| in the 
particle-hole channel, and the finite-range Gogny interac- 
tion [35| in the pairing channel. The present calculation 



of the isoscalar monopole response is carried out using the 
density-dependent meson-exchange effective Lagrangian 
DD-ME2 [H. 



B. Results and interpretation 

The values of the nuclear incompressibility calculated 
using the CHFB, non-relativistic QRPA, and relativistic 
QRPA for the chains of Sn and Pb isotopes are displayed 
in Fig. [TJ With the GMR excitation energies i^cMR 
determined by the moments of the strength distributions 
(cf. Eq. ^), Ka is calculated from Eq. ^. We note that 
the three microscopic models predict remarkably similar 
isotopic trends, with values of ¥^a for individual nuclei 
that do not differ by more than 10%. The nuclear in- 
compressibility Ka is rather constant for stable nuclei 
and systems with moderate neutron excess, but displays 
a sharp decrease beyond the neutron shells iV = 82 for 
Sn, and N = 126 for Pb isotopes. This behavior of Ka 
was already observed in Ref. |6|, in which the isotopic 
Pb and Sn chains were calculated in the QRPA with both 
the SLy4 and UNEDFO effective interactions. The rather 
constant values of Ka shown in Fig. [1] are close to 140 
MeV, in quantitative agreement with the HO prediction 
(Eq. PT7)) V This shows that the analytical approach of 
section|TT]B.2 provides an accurate approximation for Ka 
in isotopes close to stability, for which valence neutrons 
typically occupy the same major shell. In the model of 
Sec. ini the energy of the monopole is taken as 2hjj, in 
accordance with the experimental trend {huj w 41 A~^/^ 
and EisGMR ~ 80 A~^/^). Many energy functional are 
characterized by values of the effective mass m* /m that 
are smaller than 1. In this case it is the RPA attractive 
residual interaction that lowers the monopole energy to 
its empirical value, even when the unperturbed energy is 
larger than 2fta;, that is, (2hujm/m*). Therefore, one un- 
derstands that the HO estimate of Ka ~140 MeV is close 
to the results obtained EDF calculations. Even though 
the latter include effects of the effective mass and the 
residual interaction, there is an approximate cancellation 
of these effects. 

The sudden decrease of K^ beyond the doubly magic 
nuclei ^'^^Sn and ^"®Pb is, of course, related to the open- 
ing of new major shells in which monopole particle-hole 
configurations are formed. This leads to the appearance 
of a soft monopole mode that increases the monopole 
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FIG. 1. (Color online) Nuclear incompressibility of the Sn and 
Pb isotopic chains, calculated using the microscopic Skyrme- 
CHFB method (circles), the Skyrme-QRPA (triangles), and 
the relativistic QRPA (squares). 
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FIG. 2. Isoscalar monopole response in Sn isotopes, calcu- 
lated using the relativistic QRPA with the DD-ME2 func- 
tional (solid), and the QRPA with the functional SLy5 
(dashed) . The response is folded with a Lorentzian of 1 MeV 
width. 



moment m_i and, consequently, decreases the deduced 
value of K^. Soft monopole modes have been predicted in 
several studies |19l - [2ll |37| . In Ref. [20] the effect of a soft 
monopole mode on the excitation energy of the GMR was 
analyzed for Ni and Ca isotopes using the Skyrme-QRPA 
model, and similar conclusions were reached. Figs. [2]and 
13] display the isoscalar monopole response in a number 
of Sn and Pb isotopes, calculated using both the non- 
relativistic and relativistic QRPA with the functionals 
SLy5 and DD-ME2, respectively. The emergence of a 
soft monopole mode for exotic nuclei far from stability is 
clearly visible, and can be related to the decrease of the 
nuclear incompressibility Ka in neutron-rich nuclei. 

The steep decrease of Ka shown in Fig. [1] raises the 
question of the effect of magicity on the nuclear incom- 
pressibility. Despite a specific prediction jll| , no evidence 
of the mutual magicity enhancement (MEM) effect for 
the incompressibility was found experimentally [13| . The 
absence of the MEM effect was further confirmed by ad- 
ditional studies [y,[i3l- However, the present results show 
that just beyond major-shell closure the nuclear incom- 
pressibility decreases because of the appearance of the 
soft monopole mode. One could, therefore, conclude that 
a magicity effect on the incompressibility is predicted, al- 
though it is not the MEM effect. 

Figs. [5] and |3] clearly exhibit the emergence of a soft 
monopole mode in the QRPA monopole response func- 
tions of neutron-rich nuclei. It might, therefore, be useful 
to calculate the values of the incompressibility Ka sep- 
arately in the low-energy region of the soft mode, and 
in the energy region of the GMR. The results are shown 
in Figs. H] and [5] where the isotopic dependence of this 
values, calculated using the QRPA with the functional 
SLy5, is plotted for a series of Sn and Pb nuclei, respec- 
tively. The boundary between the low-energy and GMR 
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FIG. 3. Same as described in the caption to Fig. [2] but for 
the Pb isotopes. 



regions is set arbitrarily at 13.5 MeV for Sn isotopes, 
and at 11 MeV for Pb isotopes. Ksoft {Kgmr) is de- 
termined from the monopole response below (above) this 
threshold. The isotopic dependence of K^oft and Kqmr 
reflects the evolution of the soft monopole mode and the 
GMR. While in the case of the Sn isotopes no definite 
trend is found, for the Pb nuclei both Ksoft and Kgmh 
show a steady decrease with increasing neutron number. 



IV. CONCLUSIONS 

The incompressibility of finite fermionic systems has 
been studied using several analytical and numerical ap- 




no 114 118 122 126 130 134 138 142 146 



FIG. 4. Nuclear incompressibility in the Sn isotopic chain cal- 
culated using the microscopic Skyrme-QRPA with the func- 
tional SLy5. The total nuclear incompressibility Ka is shown 
in the upper panel, whereas the values determined by the 
monopole response in the low-energy and GMR regions are 
plotted in the middle and lower panels, respectively. 
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proaches. Analytic relations obtained using a simple free 
Fermi gas or harmonic oscillator (HO) models point to 
a direct link between the incompressibility of the system 
and the zero point kinetic energy of constituent particles. 
The incompressibility is thus related to the localization of 
the individual fermions. For atomic nuclei the HO model 
predicts the incompressibility Ka ~ 5Tq « 140 MeV. 

To quantitatively analyze the evolution of nuclear in- 
compressibility in isotopic chains of medium-heavy and 
heavy nuclei, microscopic EDF-based approaches have 
been used: the constrained HFB method and the QRPA 
based on both non-relativistic and relativistic density 
functionals. All microscopic models predict a remark- 
ably similar trend for the incompressibility: an almost 
constant behavior in stable nuclei and systems not far 
from stability, with values of Ka close to those predicted 
by the HO model, and a steep decrease of Ka in nuclei 
with pronounced neutron excess, caused by the emer- 
gence of a soft monopole mode in neutron-rich nuclei. 
It would be important to have an experimental confir- 
mation of the soft monopole mode in neutron-rich nu- 
clei and, more generally, to study the behavior of the 
monopole response in exotic nuclei and/or beyond shell 
closures. For instance, measurement of the GMR beyond 
^■^^Sn and ^'^^Pb would be very interesting. 
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